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The Crew Exploration Vehicle Parachute Assembly System (CPAS) is being designed to 
land the Orion Crew Module (CM) at a safe rate of descent at splashdown. Flight test 
performance must be measured to a high degree of accuracy to ensure this requirement is 
met with the most efficient design possible. Although the design includes three CPAS Main 
parachutes, the requirement is that the system must not exceed 33 ft/s under two Main 
parachutes, should one of the Main parachutes fail. Therefore, several tests were conducted 
with clusters of two Mains. All of the steady-state rate of descent data are normalized to 
standard sea level conditions and checked against the limit. As the Orion design gains 
weight, the system is approaching this limit to within measurement precision. Parachute 
“breathing,” cluster interactions, and atmospheric anomalies can cause the rate of descent to 
vary widely and lead to challenges in characterizing parachute terminal performance. An 
early test had contradictory rate of descent results from optical trajectory and Differential 
Global Positioning Systems (DGPS). A thorough analysis of the data sources and error 
propagation was conducted to determine the uncertainty in the trajectory. It was discovered 
that the Time Space Position Information (TSPI) from the optical tracking provided 
accurate position data. However, the velocity from TPSI must be computed via numerical 
differentiation, which is prone to large error. DGPS obtains position through pseudo-range 
calculations from multiple satellites and velocity through Doppler shift of the carrier 
frequency. Because the velocity from DGPS is a direct measurement, it is more accurate 
than TSPI velocity. To remedy the situation, a commercial off-the-shelf product that 
combines GPS and an Inertial Measurement Unit (IMU) was purchased to significantly 
improve rate of descent measurements. This had the added benefit of solving GPS dropouts 
during aircraft extraction. Statistical probability distributions for CPAS Main parachute 
rate of descent and drag coefficient were computed and plotted. Using test data, a terminal 
rate of descent at splashdown can be estimated as a function of canopy loading. 
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I. Introduction 

T HE Crew Exploration Vehicle Parachute Assembly System (CPAS) is required to safely land the Orion Crew 
Module (CM) at a rate of descent not to exceed 33 feet per second at sea level on a standard day. A series of 
flight tests are being conducted to test the design against this requirement. A Generation (Gen) I test called CDT-3 
appeared to possibly break this limit, depending on which data source was used. Any additional parachute material 
would add mass to Orion, and cascade into the entire launch architecture at considerable cost. A thorough analysis 
of the accuracy of each data source was therefore warranted to determine whether the requirement was met, 
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especially as the baseline Orion design gained mass. The program strives to minimize CPAS mass necessary to meet 
requirements. 

All tests discussed here were conducted at the Robby Drop Zone (DZ) at Yuma Proving Ground (YPG) in 
Yuma, AZ,* except for the Pad Abort One test at White Sands Missile Range (WSMR) in New Mexico. Each test 
was a cluster of two or three CPAS Main quarter spherical ringsail parachutes with a reference diameter, D 0 , of 1 1 6 
ft. By normalizing Main steady-state data to a standard day at sea level, CPAS can take advantage of all the valid 
data collected for similar tests. Given enough time, a functioning parachute system will “forget” its initial state, and 
the system can be defined as ergodic. Under this assumption, the velocity at test vehicle touchdown has no more 
predictive value than any other relevant instant. In fact, the tests were conducted in a desert environment with 
different ground effects than the ocean. 

The steady-state drag coefficient and rate of descent vary considerably during the Main parachute phase due to 
parachute “breathing,” cluster interaction, and atmospheric anomalies such as updrafts and wind gusts. Therefore, 
understanding parachute performance involves a combination of improved fidelity of modeling and a large base of 
tests from which to draw statistics. 


II. Data Reduction Method 


A flow diagram of the typical test 
data process is shown in Figure 1. 
Information is first assembled into 
“Best Estimate” Atmosphere, Wind, 
and Trajectory files. This paper will 
focus on the steady- state performance 
output parameters. 

The data reduction process also 
propagates instrumentation 

uncertainty through to calculated 
parameters. Each resulting function, y, 
from n uncorrelated measured 
variables, x i? can be organized in the 
form of Eq. (1). The uncertainty 
estimate, Ay, is then calculated 
according to Eq. (2), where the 
uncertainty in each measurement, Ax i? 
is multiplied by the partial derivative 
of the function with respect to that 
variable. The root sum of the squares 



Figure 1. Data Post-Processing Overview. 


is taken for each variable to compute the propagated uncertainty. (1) For consistency, all uncertainties are calculated 
at a 95% probability. 


N 


y = f(x 1 ,x 2 ,...x n ) 


1 1 / 2 


trl 


a) 

( 2 ) 


A. Position and Inertial Velocity 

The vertical position of the test article at any time will need to be converted to Mean Sea Level (MSL) to relate 
to atmospheric data from external sources. Aerodynamic vehicle and parachute performance are defined in terms of 
airspeed relative to the air mass. However, most velocity measurements are made relative to an inertial frame. Local 
winds will therefore have to be estimated and subtracted out from inertial velocity. 


* The DZ is referred to as Robby when the aircraft flies Northward, as was the case in all these tests, or as La Posa if 
the aircraft flies Southward. 
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Orthometric 



Figure 2. Conversion from Height Above Ellipsoid to Mean Sea Level. 


1. Global Positioning System (GPS) 

GPS altitude measurements are 
usually recorded as a Height Above 
Ellipsoid (HAE), h. In order to convert 
to an orthometric height or Mean Sea 
Level (MSL), H, the local geoid 
height, N, must be subtracted from the 
HAE. This equation is illustrated in 
Figure 2. (2) The geoid data were 
interpolated from tables at every point 
using databases and software obtained 
from NOAA. (3) Throughout the 
conterminous U.S., the geoid lies 
below the WGS 84 ellipsoid, so the 

geoid height is negative. (4) The geoid height at YPG is approximately -33 meters or -100 feet but can vary by several 
meters based on the exact location. 

The uncertainty of a GPS measurements are complex, which can be simplified by making some assumptions 
based on the system description for a typical test. GPS position is determined by “pseudorange” computations based 
on the time a signal travels and knowledge of the satellite position from orbital parameters. The uncertainty at a 
given instant is the product of the standard deviation of the pseudorange measurement error, cr, and a dimensionless 
term which depends on the satellite geometry called Dilution of Precision (DOP). The position error of each GPS 
device was reported by each manufacturer as a Circular Error Probable (CEP), which is defined as the radius of a 
horizontal circle which would contain the true position with a 50% probability. CEP uncertainty can be converted to 
a 95% probability by multiplying times 2.1. (5) Optimum satellite positioning may have a DOP of 1, though values of 
2 or 3 would be more common. DOP values of 5 or 6 are generally considered unacceptable for a GPS unit. The 
Vertical Dilution Of Precision (VDOP) is always higher than Horizontal Dilution Of Precision (HDOP) because the 
Earth blocks satellites that would be ideal for a vertical pseudorange estimate. VDOP can be assumed to be roughly 
twice HDOP. (4) In the absence of instant estimates, HDOP will be assumed to be 2 and VDOP will be assumed to be 
4 for a GPS with satellite lock. The resulting horizontal position uncertainties are presented in Eq. (3) and Eq. (4), 
and altitude uncertainty is given in Eq. (5). 


|AX North | = 2.1 -CEP -HDOP 

(3) 

AX East = 2 - 1CEP ' HD0P 

(4) 

|AH| = 2.1 CEP VDOP 

(5) 


The YPG facility includes GPS ground stations to improve the position accuracy. The instantaneous 
measurements at known positions are used to calculate and subtract out the inherent errors of mobile GPS receivers 
using double differencing. This method is known as Differential GPS (DGPS). Any GPS receiver can act as ground 
station provided it is left stationary long enough to obtain a highly accurate position fix and that it be recording 
simultaneously with the mobile receivers. (6) It is also possible to use the Continuously Operating Reference System 
(CORS) to improve the position solution. However, the nearest CORS station to the Robby DZ is about 40 km 
away, which is just on the edge of usability. 

GPS may determine velocity by two different methods. The relative velocity between a satellite and receiver can 
be computed from observed Doppler shift of the signal frequency. Knowledge of the satellite orbit can then be used 
to compute inertial velocity. The velocity uncertainty for most GPS models using this method is reported as an RMS 
value which is equivalent to a single standard deviation probability of about 68%. Multiplying the RMS value by 1.7 
will convert it to a 95% probability. (5) The horizontal components of uncertainty are shown in Eq. (6) and Eq. (7). 
Although documentation is lacking on this subject, it is expected that vertical velocity uncertainty will be twice as 
high as horizontal uncertainty for the similar reasons as position uncertainty. The vertical velocity uncertainty 
calculation is given in Eq. (8). (4) 
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|AV North | - 1 .7 • |AV rms 
|AVEas t | = l-7-|AV RMS | 
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( 6 ) 

( 7 ) 

( 8 ) 


CPAS used redundant NovAtel Superstar II DGPS units for the test article in Generation I testing. These provide 
inertial position and velocity (using the Doppler shift method) at rates up to 5 Hz. (7) Similar units are still used to 
track the parent aircraft, chase helicopter, and any Windpacks. 

A second method is to use an Inertial Measurement Unit (IMU) to integrate accelerometer readings once for 
velocity and twice for position. An IMU is particularly useful in airborne parachute testing where the GPS solution 
is typically lost for about 30 seconds after the test vehicle is extracted from the parent aircraft. (8) This happens as the 
GPS transitions from receiving a re-radiated signal inside the cabin to direct satellite observations. IMU errors will 
accumulate with time due to bias and drift. The optimal method is to regularly correct IMU propagation with GPS 
measurements. (9) An IMU that runs without GPS updates may not provide a sufficient level of accuracy. (10) An 
integrated GPS/IMU may provide position and velocity uncertainty estimates based upon the residuals of a Kalman 
filter. 

Two models of integrated GPS/IMU were used by CPAS beginning with Generation II testing. The NovAtel 
SPAN-SE (Synchronized Position Attitude & Navigation) is generally considered as “truth” for test vehicle position, 
velocity, and attitude. A Kalman filter in the post-processing software is used to continue position and velocity 
estimation during any GPS dropout/ 1 The attitude sensors use laser ring gyros which are not susceptible to 
magnetic interference of the metallic aircraft and payload, which is a common problem with other attitude 
sensors. (12) Due to high unit cost, tests only use a single SPAN-SE on the test vehicle. Therefore, CPAS also began 
using several Crossbow NAV440 navigation systems as backup. Each NAV440 provides full inertial attitudes, 
angular rates, accelerations and GPS position. The velocity is “inertially derived” by integrating accelerometer 
readings. (13) As a consequence, the velocity uncertainty is several times higher than GPS units that measure direct 
frequency phase shift measurements. 
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Figure 3. Sample azimuth 
DZ KTM location. 


CfbpfcU jffc-Ra** 
readings from Robby 


2. Time Space Position Information (TSPI)Optical Solution 
YPG has a series of fixed ground cameras on Kineto 
Tracking Mounts (KTM) which follow the test article or 
parent aircraft trajectory. The YPG Data Computation Team 
analyzes video from KTM cameras to determine a Time Space 
Position Information (TSPI) of payload and parent aircraft as a 
function of time. These are calculated from azimuth and 
elevations histories from each KTM boresight. The relative 
positions of the test vehicle center of mass are tracked 
manually in each video frame for precise azimuth and 
elevation fixes. A minimum of two perspectives are needed to 
compute the distance to the target at the intersection (or 
closest approach) of two rays. In practice, at least three KTM 
cameras are used to reduce error through a least squares 
method. Figure 3 shows the KTM locations at Robby Drop 
Zone (DZ). Four sample line-of-sight vectors between each 
KTM and target illustrate the azimuth readings, which are 
measured clockwise from North. 

The TSPI heights and KTM elevation angles are illustrated 
in Figure 4. The best fit TSPI location of the target is 
calculated in a Cartesian frame relative to a test origin which may change from test to test. Knowledge of the local 
datum is needed to convert the TSPI vertical coordinate to a height above mean sea level. 

YPG takes care to minimize errors in the TSPI solution by 
performing regular surveys of the KTM locations and 
measuring the site elevation and azimuth bias contributions 
before each flight test. The reported TSPI solution is smoothed 
to reduce noise. The unsmoothed minimum error coordinates at 
each time are recorded in the TSPI “Odle solution.” This also 
records residuals in azimuth, elevation, and distance from each 
KTM. The residuals can be used to estimate position 
uncertainty through extensive trigonometry. The uncertainty in 
each reading for a given KTM will sweep out a wedge-shaped volume, which can be converted to Cartesian 
coordinates. The uncertainty readings from each KTM used on test CDT-1 are shown for a single point along the 

trajectory in Figure 5. The uncertainty volumes for some KTMs are larger than others. The smoothed TSPI 

trajectory solution is nearly identical to the Odle solution, while both are bracketed by two DGPS units on the load. 

TSPI is generally an accurate measurement 
of test vehicle position and was the only 
available data source during each GPS dropout 
for Gen I tests. The first Gen II test to solve the 
dropout problem by using a GPS/IMU was 

called EDU-A-TSE-1A. The altitude 

measurements from each source for this test are 
compared in Figure 7 with vertical uncertainty 
estimates shown close up. The TSPI uncertainty 
is calculated as an average of the vertical 
uncertainties from each KTM, and can be seen 
to be more accurate than DGPS measurements. 
The SPAN-SE data is considered the best 
available data though the indicated error bars 
are overly conservative estimates based on the 
specification sheet, because the Kalman filter 
error estimates were not available for this test. 

Figure 5. TSPI uncertainty volumes for CDT-1. There is also some small discrepancy in the 

datums used in post-processing each system. 

The uncertainty bounds overlap, implying that a true physical solution is contained between them. 



Figure 4. Sample KTM elevation readings. 
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Despite its precision with position, TSPI is not 
the optimum method to compute velocity, because it 
must be determined indirectly through numerical differentiation. The rate of descent from TSPI, for example, is the 
time rate of change of altitude. To reduce errors, a smoothed solution is taken through the nominal data. In the same 
way, two other solutions can be made by calculating the slopes through the minimum and maximum altitudes, as 
illustrated in Figure 7. The difference between these velocities can be used as a velocity uncertainty estimate at each 
point, using Eq. (9). 


Figure 6. GPS and TSPI altitude for EDU-A-TSE- 1 A. 
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A sample of the downward velocity and error 
estimates from multiple data sources is shown in 
Figure 8 for test EDU-A-TSE-1A. The GPS/IMU 
tends to track well with both DGPS units. The 
uncertainty estimates for TSPI velocity are much 
larger than for the GPS systems. Further, the GPS 
systems tend to show a more damped trajectory, 
which is more realistic based on the physics of the 
situation. As with the earlier Gen I two-Main 
cluster test, the rate of descent from this test comes 
very close to the design limit, justifying the 
acquisition of the higher accuracy system to 
optimize the design. 


EDU-A-TSE-1A, October 2, 2009 



Figure 8. Vertical velocity comparison for EDU-A-TSE- 1 A. 
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B. Atmospheric Measurements 

In order to compare Main parachute steady-state performance across the entire altitude range within a test or 
from one test to another, the data are normalized to sea level equivalent conditions based on ambient atmospheric 
measurements. These measurements are provided by Radar Wind Sounding (RAWIN) balloons and a ground 
weather station. The uncertainty of the atmospheric measurements is combined with the uncertainty of the inertial 
velocity measurements through a series of uncertainty propagation equations. Temperature, T, static pressure, p, and 
relative humidity, RH, are directly measured to calculate ambient air density, p. The saturation pressure, p sat , is first 
determined with Eq. (10) where the temperature must be in Kelvin and the result is in mb. The uncertainty in 
saturation pressure is given by Eq. (11). The humidity-corrected density is computed in Eq. (12). Using the four 
independent variables p, RH, p sat , and T, the uncertainty is computed by taking the partial derivative of the density 
equation with respect to each, according to Eq. (2), which simplifies to Eq. (13). 


p sat =6.1078-10 


7-5-1V. 


-2048.625 


Kelvin “35.85 


( 10 ) 


| A Psat | = |p sat ( T Kelvin + | AT |)- Psat ( T Kelvin " | AT || 


(11) 


P • p sa t ^ RH • p sat 

R T R T (12) 

AV air A AV vapor A v y 


\Ap\ = 

r n ^ 

2 

_i_ 

( 

ARH • 

R, T , 

i 


air J 


V v 


P sat _j_ P sat 




R . T R • T 

^air ’ 1 ^ vapor 1 J J 


| A Ps 


f -RH RH 1 

i 

\2 

+ 

/ 

FT 

< 

2 

1/2 

i 

R • T R . t 

^^air ' 1 ^ vapor ’ 1 j 

l T J 


(13) 


C. Steady-State Parachute Performance 

A single parachute is considered to be in steady-state descent when 
the drag force is equal to the payload weight. The force balance for a 
cluster is more complex, as illustrated in Figure 9. Steady-state can be 
determined when the test vehicle has first decelerated to 1 g. 

The primary source of the wind profiles are Windpacks. These are 
instrumented payloads suspended under very stable tri-lobe parachute 
system designed for a predictable equivalent rate of descent. Any 
horizontal displacement is attributed to the local wind field. (14) Because 
the DGPS units in Windpacks experience a dropout at aircraft release, the 
RAWIN balloon data are used to fill in this gap. Ground station readings 
are also included in the Best Estimate Wind profile. Horizontal 
components of payload airspeed are computed by subtracting the 
horizontal wind components from the inertial velocity. Instantaneous 
updrafts and downdrafts are currently not measured due to their transient 
nature. The inertial vertical velocity and wind-corrected horizontal 
velocity are used to compute local airspeed, V air . Riser tension can be 
decomposed by measuring parachute fly-out angles relative to the 
geometric centroid. This allows for comparison with loads from an 
accelerometer. 

Simplifications are used to compute the cluster steady-state drag 
coefficient in Eq. (14). The total parachute area is the number of 

Figure 9. Parachute cluster diagram. 


Parachute Drag & Normal F orces 


Centroid^.v M 


Fly-out \ 
Angles, 0 X \ 



q-(C D S) 


Payload Weight, 
W 


Vertical 
Velocity, 


Horizontal 
Wind 
Field 


V- 


American Institute of Aeronautics and Astronautics 



parachutes in the cluster, N c , multiplied by the reference diameter, S 0 . By CPAS convention, the full open steady- 
state drag coefficient, C Dq , is normalized to the suspended weight of the payload, W Pay i oad , which does not include 

the weight of the parachute material and suspension equipment. The parachute steady-state drag coefficient is 
defined using purely vertical velocity, V z , which neglects any parachute gliding motion. To isolate parachute drag, 
the drag area of the payload or vehicle, (C D S) Pay i oa d, must be subtracted from the system drag area. 


Cn - 
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N P -S n 


^Payload 
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v2 z 


-(C D S) 


Payload 


(14) 


Equation (14) has four independent variables: W Payk)ad , V z , p, and (C D S) Pay i oad . The partial derivative of the drag 
coefficient equation is taken with respect to each of these using the standard uncertainty analysis. This leads to the 
simplified drag coefficient uncertainty in Eq. (15). 
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Parachute rate of descent performance at altitude with varying atmospheric conditions are normalized to the sea 
level equivalent rate of descent, V ZsL , for a standard day as defined in Eq. (16). 


Vy 



(16) 


Uncertainty in equivalent rate of descent is computed according to Eq. (17). Uncertainty for both drag coefficient 
and equivalent rate of descent are most sensitive to vertical velocity, which is one of the reasons CPAS has invested 
in high quality GPS instrumentation. They are also both sensitive to density, which is not measured directly. 
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The most important driver of rate of descent is the canopy loading, defined as the suspended weight divided by 
the total parachute surface area, W pay i oad /(N c -S 0 ). A summary of test results vs. canopy loading is presented in the 
conclusion. 

Note that both C Dq and V Z§L parameters are based on the same test data. Drag coefficient is proportional to the 

inverse square of equivalent rate of descent, so for a given test condition, a curve of approximately Y=l/X 2 can be 
generated. Equations (14) and (16) can be combined to make the relation in Eq. (18). 
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CPAS parachute test data suggest that the distribution of sea level equivalent rate of descent data is Gaussian. If 
that is the case, then it follows from Eq. (18) that the corresponding distribution for drag coefficient is skewed with a 
peak probability density lower than the mean and a long tail towards high drag (or low rate of descent) end of the 
curve. Care must be taken when referring to statistical terms such as mean and standard deviation, which may not 
apply to drag coefficient. 
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III. Selected Flight Test Results 

Flight data are presented for selected CPAS tests, as well as the first Pad Abort test. Only tests with a cluster of 
Main parachutes using the baseline design are considered. CPAS Gen I and Gen II tests are summarized in Ref. (15) 
and Ref. (16), respectively. The results of additional design changes are discussed in separate papers. 


CDT-1, October 18, 2007 


£ 40 
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> 35 


Best Estimate Trajectory 
— “ “ Limit for Two Mains: 33 ft/s 
Time Average: 21.34 ft/s 


A. CDT-1 

CDT-1 was conducted on October 18, 2007. The test 
vehicle was a parachute compartment mockup mounted on a 
9x20 ft Type V platform. Extraction from the C-130A occurred 
at about 10,930 ft MSL where the system immediately 
deployed two CPAS Drogue parachutes. Pilot parachutes were 
mortar deployed to successfully extract the cluster of three 
Main parachutes already shown in Figure 9. Steady-state 
descent began at about 3,790 ft MSL and lasted for about 110 
seconds. The suspended weight was set to approximate the 
Orion design at the time of 16,462 lb. With three Mains, the 
canopy loading was 0.519 lb/ft 2 . 

A time history of the sea level equivalent rate of descent is 
shown in Figure 10. Due to the fairly low canopy loading, the 
system never approached the performance limit. The drag 
coefficient shown in Figure 11 showed an increasing trend. 

This is an artifact of reducing the data as if the airspeed was a 
purely vertical, when in fact the parachutes were gliding. 

Further, a temperature inversion layer was encountered. This is due to the standard procedure of testing soon after 
dawn when winds are minimal, yet the desert floor begins warming. 

A statistical analysis of these data is presented in Figure 
12. A histogram of equivalent rate of descent is shown in the 
upper right. The data were fit as a Gaussian curve, shown in 
black, to calculate a time average of 21.34 ft/s and standard 
deviation, cr, of 1.25 ft/s. Bounds at the mean plus and minus 
3cr are indicated, which slightly exceed the minimum and 
maximum data actually encountered. These points are 
mapped to the drag coefficient histogram through the inverse 


Time after Ramp Clear, 14:42:50.47 (s) 


Figure 10. CDT-1 equivalent rate of descent. 


CDT-1, October 18, 2007 


Time after Ramp Clear, 14:42:50.47 (s) 

Figure 11. CDT-1 drag coefficient. 

bounds are not described as 3 cr, but as 
the equivalent probability of 99.73%. 

Randomized rate of descent data 
were generated using the calculated 
mean and standard deviation. This curve 
is plotted in red. These random data 
were then mapped to the drag 
coefficient histogram on the lower left 
and plotted as a red curve. This curve far 
better matches the drag coefficient 
distribution than attempting to fit drag 


square relationship plotted in the lower right. If the 
rate of descent CDT-i,ocoteM8,2007 

distribution is truly 
Gaussian, then the 
drag coefficient 
cannot 

simultaneously be 
Gaussian. Therefore, 
the drag coefficient 





coefficient 

distribution. 


with 


a Gaussian Figure 12. CDT-1 steady-state performance histograms. 
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The drag coefficient can be described as “right skewed,” that is the highest probability density is lower than the 
time average and a long tail points toward high drag coefficient. The time average of drag coefficient is 0.962, but 
the time average rate of rate of descent corresponds to 0.951. Designing the parachute using the lower drag 
coefficient is the more conservative approach. 


B. CDT-3 

CDT-3 was conducted on June 17, 2008 and employed a cluster of two CPAS Mains 
seen in Figure 13. CPAS is required to meet the landing requirement with as few as two 
Mains, so several tests were conducted to drive the design with this limiting case. The test 
vehicle was a weight tub mounted on a 9x20 ft Type V platform. Extraction from the C- 
130A occurred at about 18,100 ft MSL where the system immediately deployed two 
CPAS Drogue parachutes. The Mains were deployed by static line as the Drogues cut 
away. 


CDT-3, June 17, 2008 




Figure 14. CDT-3 equivalent rate of descent. 


Steady-state descent began at 
about 6,400 ft MSL or about 130 
seconds before touchdown. The 
suspended weight was about the 
same as CDT-1 at 16,717 lb but the 
canopy loading was higher at 0.719 
lb/ft 2 due to using one less Main. A 
strong temperature inversion layer 
was encountered at about 3,200 ft 
MSL, and vertical winds were Figure 13. CDT-3 
indicated by both Windpacks until the surface. The assumption 
of normalizing drag coefficient by the inertial vertical velocity 
breaks down when significant gliding or vertical winds are 
present. Therefore, statistics were only taken until about 223 
seconds after ramp clear such that the time average rate of 
descent is 26.73 ft/s. 

The equivalent rate of descent is plotted in Figure 14. The 
nominal TSPI data exceeds the 33 ft/s limit on two occasions. 
The second of these occasions is shown in more detail for each 
of the raw data sources. The uncertainty of TSPI is 
consistently much higher than the uncertainty from DGPS 
readings. For most of the trajectory, there is overlap of these 
uncertainties, which validates the TSPI uncertainty estimate 
method and DGPS specification sheets. Therefore, the Best 
Estimate Trajectory (BET) is constructed from the average of 
the two DGPS units, neglecting the TSPI. This investigation 
indicates that the payload likely never actually broke the limit, 
and that the parachute design was successful. The 
corresponding Main steady-state drag coefficient is plotted in 
Figure 15. 
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CDT-3, June 17, 2008 
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Figure 15. CDT-3 drag coefficient. 


Histograms of steady-state performance are shown in 
Figure 16. Discounting the terminal data past encountering 
the inversion layer makes the histograms more Gaussian. 
The fitted standard deviation to rate of descent is 1.22 ft/s 
and the weighted drag coefficient is 0.921. 


CDT-3, June 17, 2008 
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Figure 16. CDT-3 steady-state performance histograms 
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C. EDU-A-TSE-1A 


CPAS Generation II testing began 
with a series of flights named EDU- 
A-TSE-1 or “smart release.” The 
primary objectives dealt with avionics 
development but these tests were a 
good opportunity to get two-Main 
cluster performance as a secondary 
objective. EDU-A-TSE-1A was 
conducted on October 2, 2009. The 
smart release vehicle was a modified 
weight tub attached to a 9x24 ft Type 
V platform. Extraction from the C- 
130A occurred at about 21,574 ft 
MSL. A programmer parachute was 
Figure 17. EDU- deployed via a delayed load transfer 
-T E- Mains. to set up the test condition. A cluster 
of two Mains, seen in Figure 17, were static line deployed as 
the programmer cut away. Black dye markings were added to 
the skirt in an attempt to improve visibility for 
photogrammetrics. Each Main parachute also incorporated an 
Over-Inflation Control Line (OICL) which was intended to 
reduce the variation of rate of descent by limiting the inlet area. 

The test article weight was chosen as 21,574 lb to 
approximate the design Orion mass. The increased canopy 
loading of 0.909 lb/ft2 led to a higher time average rate of 
descent of 28.49 ft/s seen in Figure 18. 

. On three occasions the nominal rate of descent from the 
NovAtel SPAN-SE approaches the 33 ft/s limit with 
uncertainty bands exceeding the limit twice. These sudden 
increases of descent rate, or sudden drops in drag coefficient in 



EDU-A-TSE-1A, October 2, 2009 



Figure 18. EDU-A-TSE-1A equivalent rate of 
descent. 


EDU-A-TSE-1A, October 2, 2009 


Figure 19, are due to relative instability of the cluster 
formation causing interference between the Mains. Like most 
of CDT-1, the drag coefficient for this test shows a slight 
downward trend. 

The OICL had negligible effect on the amount of 
variation of rate of descent. The standard deviation of 
equivalent rate of descent for this test was 1.50 ft/s, which 
was higher than the CDT-3 value of 1.44 ft/s. A 
photogrammetric analysis showed that the skirt perimeter 
may never have completely opened to the OICL length. 

A statistical analysis of these data is presented in Figure 
20. A histogram of equivalent rate of descent is shown in the 
upper right. The standard deviation of the fitted rate of 
descent is 1.50 ft/s. The flight test data does not actually 
encounter the mean plus and minus 3 <j equivalent values. 
The randomized fit curves in red generally describe histograms. The typical right skew in the drag coefficient is not 
seen in this case as there were no exceptionally high drag data points. 



Figure 19. EDU-A-TSE-1A drag coefficient. 
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Drag Coefficient, C 


EDU-A-TSE-1A, October 2, 2009 




Figure 20. EDU-A-TSE-1A steady-state performance histograms. 
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D. EDU-A-TSE-1B 



EDU-A-TSE-1B, December 1, 2009 


r 


Figure 21. 
1A Mains. 


EDU-A-TSE- 


Best Est. Traj. (NovAtel GPS/IMU) 
— — - Limit for Two Mains: 33 ft/s 
Time Average: 28.56 ft/s 


28 ^ - #1 -f W§ 


Time after Ramp Clear, 1 5:21 :38.4889 (s) 


Figure 22. 
descent. 


EDU-A-TSE-1B equivalent rate of 


The second smart release 
test, EDU-A-TSE-1B was 
conducted on December 1, 

2009. The same test vehicle 
and concept of operations was 
used as the previous test. 

Extraction from the C-130A 
occurred at about 21,175 ft 
MSL. The cluster of two mains 
successfully deployed as seen 
in Figure 21. Additional black 
markings were added to 
different colored panels near 
the crown to aid in parachute 
identification of upward-looking video. 

A time history of the sea level equivalent rate of descent is 
shown in Figure 22. The drag coefficient is shown in Figure 
23. The high descent rate and low drag event near the 95 s mark was due to the parachutes flying out and crashing 
against each other, lowering the effective canopy surface area. 

A statistical analysis of these data is presented in Figure 24. A histogram of equivalent rate of descent is shown 

in the upper right. The data were fit as a Gaussian curve, 
shown in black, to calculate a time average of 28.56 ft/s and 
cr of 1.44 ft/s. The test data never went past the mean V Z§L - 

3 <t bound, but do exceed the mean V Zsl + 3 cr bound due to 

the previously described event. The randomized red curves 
generally fit the shape of 
the histogram data with the 
exception of tail created by 
that region. 

The time average of 
drag coefficient is 0.931, 
but the time average rate 
of rate of descent 


EDU-A-TSE-1 B, December 1 , 2009 
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Figure 23. EDU-A-TSE-1B drag coefficient. 

corresponds to a more conservative drag 
coefficient of 0.924. 
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Figure 24. EDU-A-TSE-1B steady-state performance. 
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E. Pad Abort One 


The first Pad Abort demonstration test, PA-1, took place on May 6, 
2010 at the White Sands Missile Range. A representative capsule shape 
ascended using the Launch Abort System (LAS) and descended under a 
sequence of CPAS parachutes, shown in Figure 25. The Generation I 
design meant the canopy loading under three Mains was only 0.5135 lb/ft 2 , 
lower than the Generation II cluster tests. Due to the relatively low altitude 
at parachute deployment, the time suspended under the Main parachutes 
was only about 74 seconds. 

CPAS was not involved in the data collection for this test, but was 
provided a BET and atmosphere data. The prime source of navigation data 
from PA-1 came from sensed accelerations and angular rates in the Space 
Integrated GPS/INS (SIGI), provided by Honeywell/ 17) Other radar and 
optical data sources were combined with the SIGI data in an External 
Kalman Filter (EKF) to generate a BET. However, the system did not have 
a GPS receiver to correct inertial drift. 

Due to its similarity with CDT-1, the steady-state results are compared 
to PA-1 in Figure 26 and Figure 27. Although the trends differ from CDT- 
. 1, the time average rate of descent is very similar at 21.69 ft/s. The shorter 

igure . - am parac utes. time for steady-state statistics may have caused the data to be less Gaussian 

than other tests, as shown in the histograms of Figure 28. 



CDT-1 compared with PA-1 
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Figure 26. EDU-A-TSE-1B equivalent rate of 
descent. 


Figure 27. CDT-1 and PA-1 drag coefficient. 


16 

American Institute of Aeronautics and Astronautics 



Drag Coefficient, C 


0.5 r 
~ 0.4 h 


PA-1, May 6,2010 


Time Avg 
^ 4 . 6 £ 


Test Data 
Fitted Test Data 
Random Fit 



o 

Q 


1.8 * V 

1 

1.6 

1.4 j- 

i 

1.2 



0.6 


6 


5 



4 3 2 1 0 

Probability Density 


i. 

i. 

i. 

i. 


o. 
o. 

Sea Level Equivalent Rate of Descent, V ZSL (ft/s) 







99.7 

3% bounc 

1: 



r.43 





Weighted 







0,912 ; ' 

99.73% bOunc 

1: 


Figure 28. PA-1 steady-state performance histograms. 
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IV. Conclusion 


A summary of all time averaged CPAS 
baseline Mains cluster steady-state sea 
level equivalent rate of descent test data is 
presented as a function of canopy loading 
in Figure 28. Horizontal error bars are due 
to uncertainty of determining the test 
vehicle suspended weight. Vertical error 
bars at the data points indicate the average 
Vz SL uncertainty for the given test. The 

minimum and maximum encountered 
values of steady-state V Z§L are also 

plotted, as well as estimated 3 <j dispersion 
envelopes. A second order polynomial 
curve fit was used to fit the time averaged 
data as a function of canopy loading. 

It is obvious from the data that the 
possible rate of descent may vary 
considerably at splashdown due to the 
complexity of cluster interaction and 
atmospheric anomalies. However, the 
curve fit, reproduced in Eq. (19), can be 
used as a predictive tool for the baseline design. Equivalent rate of descent appears to have a Gaussian distribution, 
so the vertical velocity at splashdown may be modeled from Monte Carlo simulations using a random Gaussian 
distribution centered on the nominal. The average standard deviation (la) from all four cluster tests of 1.40 ft/s is 
recommended to define the distribution. Note that this equation only models the vertical velocity, so high winds may 
increase the resultant total landing velocity. 



Figure 29. 
Summary. 
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(19) 


At the heavier canopy loading, the rate of descent is likely to approach or break the limit for two Mains of the 
current design. Therefore, further CPAS tests focused on changing the Main parachute design for improved cluster 
stability. This was done by increasing the geometric porosity and lengthening the suspension lines. Further tests may 
possibly use a combination of the two. 
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